Chiral Phase Transition beyond Mean Field Approximation 



Lianyi He and Pengfei Zhuang 

Physics Department, Tsinghua University, Beijing 100084, China 
(Dated: February 1, 2008) 

Based on the analogy between the Nambu-Jona-Lasinio model of chiral symmetry breaking and 
the BCS theory of superconductivity, we investigate the effect of qq pair fluctuations on the chiral 
phase transition. We include uncondensed qq pairs at finite temperature and chemical potential 
in a self-consistent T-matrix formalism, the so-called GoG scheme. The pair fluctuations reduce 
significantly the critical temperature and make quarks massive above the critical temperature. 
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It is well known that the spontaneous breaking of 
chiral symmetry is one of the most important features 
of the vacuum of Quantum Chromodynamics (QCD). 
The spontaneous breaking successfully explains many low 
energy phenomena in QCD, such as the large contin- 
ues quark mass and the small pion mass. The quark- 
antiquark condensate (qq) which characters the sponta- 
neous breaking is about (250 MeV) 3 from the QCD sum 
rule and lattice calculations. It is generally believed that, 
in hot and dense medium, the chiral condensate will be 
reduced and the broken chiral symmetry will be restored 
at sufficiently high temperature and density. 

Since the chiral symmetry breaking and restoration 
happen in the non-perturbative region, it is hard to han- 
dle them directly by the original QCD lagrangian. The 
investigation relies mainly on low energy effective models 
and lattice calculations. One of the successful models is 
the Nambu-Jona-Lasinio (NJL) model [l|, |2| . The idea 
of the model is from the well-known Bardeen-Cooper- 
Shriffer (BCS) theory of superconductivity. In supercon- 
ducting metals, the Cooper pairing between electrons in 
the spin-singlet state leads to a condensate (V'fV'i) which 
spontaneously breaks the electromagnetic U(l) symme- 
try and causes an energy gap. In analogy, the pairing 
between a quark and an antiquark with the same chi- 
rality leads to a condensate (qq). Such a condensate 
spontaneously breaks the chiral symmetry and gives a 
dynamical mass to the quarks. 

While the analogy between the QCD and BCS vacua 
leads to a successful theory of chiral symmetry break- 
ing, one should keep in mind that there is an important 
difference between the two: In contrast to the BCS su- 
perconductivity, the chiral symmetry breaking is a strong 
coupling phenomenon. From the recent studies on rela- 
tivistic heavy ion collisions, there may exist a strongly 
coupled quark-gluon plasma phase (sQGP) [3| above the 
critical temperature T c of the deconfinement and chiral 
phase transitions where both quarks and their bound 
states are constituents of the system. From lattice sim- 
ulations, any thermodynamic quantity can not reach its 
Stefan-Boltzmann limit even at extremely high tempera- 
ture Q. A possible explanation for this thermodynamic 
suppression is that quarks have a large thermal mass in 
the deconfined and chiral restored phase [S IS @] • 

In the language of condensed matter physics, these 



phenomena strongly indicate that there exists a pseu- 
dogap for quarks [7], H, H, EH and the matter may be in 
the crossover region from BCS to Bose-Einstein conden- 
sation (BEC) p], HI EH- It is generally believed that 
the BCS mean field theory describes well the crossover 
at zero temperature. However, to obtain a quantitatively 
correct T c and a correct fermion excitation above T c , one 
should properly consider pair fluctuations at finite tem- 
perature. There are some T-matrix schemes to treat the 
finite temperature pair fluctuations. Among them, the 
self-consiste nt p air approximation or the so-called GoG 
scheme [lU, [3 EH is a competitive one. The effect of 
the pair fluctuations in this scheme is treated as a pesu- 
dogap. Recently, the GoG theory is applied to the study 
of color superconductivity in the NJL model [16( . In this 
paper, we will investigate the chiral phase transition in 
the NJL model, by using the self-consistent GoG theory. 

To have a BCS-like description of the spontaneous 
chiral symmetry breaking, it is useful to take the two- 
dimension Wely spinors defined as 



Q = 



9B 



1 = ( 9r Ql ) 



(1) 



where the color and flavor indexes of the quark field are 
not explicitly shown. The kinematic term of the NJL 
Lagrangian, Cq — qijQq, can be written as 



Co = q{i(d t - a • V)g L + qli(d t + er ■ V)g R , (2) 

where <Ji are the Pauli matrices. The interaction between 
quarks and antiquarks in the NJL model are expressed by 
a four fermion coupling term, it can be generally written 



as 



r - 9 



(qqf + (qiisFq)' 



(3) 



where Y = 1 and <r are, respectively, for one and and two 
flavor cases. To see directly the analogy between the NJL 
model of chiral symmetry breaking and the BCS theory 
of superconductivity, it is intuitive to look into the one 
flavor case. Using the identities for the Wely spinors, 



qq 
qiisq 



9r?l + qiqn 



-i(qkqt 



(4) 
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the four fermion interaction can be written as 



A = gqRqhqtqR, 



(5) 



which has the same structure of the interaction between 
fermions in BCS theory where one has C\ = gip^li'li't- 
The order parameter field describing chiral symmetry 
breaking can be defined as $(t, x) = — .gg R QL- The QCD 
vacuum is characterized by the chiral condensate ($) 
which means the pairing between a quark and an an- 
tiquark with the same chirality. Since the vacuum has 
a definite parity, the chiral condensate is a real number, 

(<7r?l) = (<7l9R>- 

We start by rewriting the mean field theory of the two 
flavor NJL model in a T-matrix formalism. Such a for- 
malism is important for us to go beyond the mean field 
by including the non-condensed qq pairs at finite temper- 
ature. The mean field quark propagator S reads 



s-\k) 



■m 



(6) 



where ^olr(^) = ^ v n + A 1 ± o - • k are the free quark 
propagators with v n being the fermion frequency v n — 
(2n + l)nT (n = 0,±1,±2, • • •), and m sc = ($) is the or- 
der parameter of the phase transition. The quark chem- 
ical potential fi is introduced by considering conserved 
charge density q^q. The explicit form of the quark prop- 
agator can be expressed as 



S(k) = 



Gh(k) T L (k) 



(7) 



with the matrix elements given by 



&,r(*) = ^olVW - s L, R (fe) 
-?L,R( fc ) = m sc Gh,R(k)GoR,-L{k), 



(8) 



where the self-energies £l,r at mean field level are de- 
fined as 



SL,R(fc) = m^ c G R,h(k). 



(9) 



After some simple algebra, the diagonal and off-diagonal 
elements can be evaluated as 



GuR(k) = 



iv n + /i =p o- • k 

(iu n + ^) 2 - k 2 - r, 



(iv n + /i) 2 



k 2 



(10) 



For those who are familiar with the BCS theory of su- 
perconductivity, G and T are in analogy to the normal 
and anomalous Green functions, and the quark chemical 
potential plays the role of effective Zeeman splitting be- 
tween quarks and antiquarks. At fi — 0, the quark prop- 
agator has two poles iv n — ±-\/k 2 + to 2 ,, which means 
that quarks obtain a mass gap m sc in the chiral symmetry 



breaking phase. The order parameter to sc is determined 
by the self-consistent gap equation 



k 



(11) 



where the summation over quark momentum is defined as 
J2k = J rf 3 k/(27r) 3 and the trace is taken in color 

and flavor spaces. 

In the mean field theory, qq pairs enter into the prob- 
lem only through the condensate. The condensed pairs 
can be associated with a T-matrix in such a way 



tsc(q) = -f^S(q) 



(12) 



that the quark self energies (|9|) can be formally expressed 



as 



El,r(A0 = ^2 *sc(9)£or,l(& - q)- 



(13) 



From the gap equation (|TT|) . the mean field theory is 
related to a particular asymmetric pair susceptibility [13| , 

x(q) = xlr(<z) = XRh(q) 

= lj2 Tl [ g ^ k ^R(k - q) 

k 

+GL(k - q)G R(k)\, (14) 

from which the gap equation in the symmetry breaking 
phase is given by the condition 



1 + 9X(0) = 0. 



(15) 



This suggests that the uncondensed qq pair propagator 
should be in the form 



* P g(<?) 



g 



i + gx(q) ' 



(16) 



The gap equation (fTS")) is just the so-called BEC condition 
tpg(q = 0) = 0. While the uncondensed qq pairs play no 
role in the BCS mean field theory, such a specific choice 
of the pair susceptibility and the BEC condition are the 
fundamental criterion for us to go beyond the mean field. 

We now take into account the uncondensed qq pairs. In 
the mean field theory, the quark self-energies El,r con- 
tain only the contribution from the condensed qq pairs. 
This is mostly correct at zero temperature. However, at 
finite temperature, the condensed pairs with zero total 
momentum can be thermally excited. Therefore, the to- 
tal contribution to the quark self-energies should include 
both the condensed (sc) pairs and the uncondensed or 
"pseudogap" -associated (pg) pairs [13j, 



£l,r(&) = ^2 i ('?)^oR,L(fc - q) 



(17) 
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with the total pair propagator t(q) defined by 



*(<?) 
t sc (q) 



tpgiq) +*sc(<?), 

9 



(18) 



With the modified self-energies, the dressed quark prop- 
agators <?L,R.(fc) and the pair susceptibility x(q) are still 
given by © and ([13]). The beyond-BCS effect is re- 
flected in the quark self-energies £l,r- The BEC con- 
dition t~g(0) — and the quark self-energies form in 
principle a coupled set of equations to determine the new 
order parameter m sc . 

The above equations are hard to handle analytically. 
In the symmetry breaking phase with T < T c , the BEC 
condition t~g(0) = implies that t pg (q) is strongly 
peaked at q — 0. This allows us to approximate [HJ 



(19) 



where m 2 includes the contribution from both the con- 
densed and the thermally excited qq pairs, 



2 2 

in = m a „ 



with the pseudogap m pg defined as 



*pg(e)< 



(20) 



(21) 



Note that, above the critical temperature T c , the BEC 
condition disappears and then such an approximation 
is in principle not valid. How to practically treat the 
pair fluctuations and their effect on the quark propaga- 
tor above T c is still an open question in the GqG scheme. 
However, for temperatures not much higher than T c , this 
approximation may be still good to give qualitatively cor- 
rect result. 

Under the approximation, the dressed quark propaga- 
tor can be analytically evaluated as 



Sl,r(*0 



H =p a ■ k 



(iv n + fi) 2 - k 2 - 



(22) 



which means that, when one goes beyond the mean field, 

) x l 2 rather 



pg 



quarks obtain a mass gap m = (m 2 c 
than to sc . One may worry about that the appearance of 
the pseudogap m pg explicitly breaks the chiral symmetry. 
Parallel to the discussion in the BCS-BEC crossover [15j. 
we can show that m pg is just the classical fluctuation of 
the order parameter field $, 



m 2 pe = (m 2 )-(m) 2 , 



(23) 



and hence does not break the chiral symmetry. 



After introducing the pesudogap m pgl the pair suscep- 
tibility can be analytically evaluated as 



rf 3 k 



1 + e— cos q 



|k - q| + eE u 



|3 (^„) 2 -(|k-q|+eSk) 2 

[l-/(eE k )-/(|k-q|)] , 
(24) 



where iV c and Nf are color and flavor numbers, and 
the pair frequency u) n , quark energy E^, the angle </> q 
and the function f(x) are defined as u> n = 2nirT (n — 
0,±l,±2,---),^k = Vk 2 + m 2 ,cos0 q = k • (k - 
q)/(|k| |k - q|) and f(x) = [f(x — /i) + f(x + fj)] /2 with 
f{x) = 1/ (e x / T + l) being the Fermi-Dirac distribution 
function. The BEC condition 1 + <?x(0) = 0, namely the 
gap equation can now be explicitly written as 



1 - gN c N f 



d 3 k l-2/(£ k 
(2tt) 3 



0. 



(25) 



The equations l[2T|) and |25|) form a closed set to deter- 
mine the order parameter m sc and the pseudogap 77i pg as 
functions of temperature and chemical potential. 

Solving such a coupled set of equations is still rather 
complicated. In the temperature region below and 
around the critical temperature T c which is expected to 
be much smaller than the momentum cutoff A of the NJL 
model, we can employ the pole approximation for the pair 
propagator t pg , 



* Pg (w,q) ^ -- 



0J„ 



(26) 



with w 2 = v 2 q 2 + 
continuation iui n - 
A 2 are defined as 



A 2 , where we have taken the analytical 
-> u!+iQ + , and the coefficients Z,v 2 and 



1 d 2 x 
1 



1 d 2 x 



u)— q— 

- 9X(0)] 



2Zdq 2 



oj— q— 



(27) 



T he p air width is safely neglected below and around 
T c [13j. With the help of the pole approximation, the 
pseudogap equation (|2"Tj) takes the simple form 



m„ = — 

Pg 7 



d 3 q l + 26(w q ) 



2uj„ 



(28) 



where b(x) — 1/ {e x l T — l) is the Bose distribution func- 
tion. Since the vacuum is well described by the mean 
field theory, we require m pg = at zero temperature and 
chemical potential. To this end, we simply eliminate the 
first term on the righthand side of ([28]). 

Before numerical calculations, we can analytically ar- 
rive at the following conclusions: 

l)In the symmetry breaking phase, the gap equation 



4 



1 + <7x(0) = leads to A 2 = 0, and then the pair dis- 
persion is gapless in the long wavelength limit. At suf- 
ficiently low temperature, the pseudogap m pg is propor- 
tional to the temperature T. 

2)From \-> % ~ N c , we have mpg ~ 1/AT C . In the large 
N c limit, the pseudogap effect can be neglected, and the 
mean field approximation describes well the chiral phase 
transition. 

The NJL model is non-renormalizable, the simplest 
regularization scheme is to use a three momentum cutoff 
A to regularize the integrals over quark and pair mo- 
menta. In the vacuum, the spontaneous breaking of 
chiral symmetry occurs when the coupling exceeds the 
critical value g c = Ait 2 / {N C N f A?) . In the following we 
use a dimensionless quantity rj — g/g c — 1 to denote 
the coupling strength. To fit the pion decay constant 
/tt = 94 MeV and constituent quark mass m = 312 MeV 
in the vacuum, we take A = 653 MeV and g/A = 5.01 
(GeV) -2 [ItJ which result in the physical coupling -q — 
0.3. In FigfTJwe calculate the transition temperature T c 
as a function of the coupling in mean field approximation 
and in the case including pseudogap effect with N c = 3 
and 100. It is clear that the pseudogap effect can be ne- 
glected at large enough N c and disappears in the limit 
N c — > oo. For a finite value of N c , the transition temper- 
ature with contribution from pair fluctuations is lower 
than the one in mean field approximation, and the differ- 
ence between the two increases with increasing coupling. 
This conclusion is consistent with the results obtained 
in other approaches, like the non- linear sigma model ap- 
proach of the NJL model 0, H|- 
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FIG. 1: The critical temperature T c of chiral phase transi- 
tion as a function of the coupling strength 77 in mean field 
approximation (dotted line) and in the case including pseu- 
dogap effect with N c = 3 (solid line) and 100 (dashed line). 
The chemical potential (x is taken to be zero. 

In the following we focus on the physical case of N c = 3 
and j] = 0.3. In Fig[2]we show the temperature depen- 
dence of the order parameter m sc , the pseudogap ra pg 
and the total quark mass m at zero chemical potential. 
At zero temperature, the pseudogap is zero and the quark 
mass comes purely from the dynamical symmetry break- 



0.5 




0.2 0.4 0.6 0.8 1 



T/T 

c 

FIG. 2: The temperature dependence of the order parameter 
m sc (dotted line), the pseudogap m pg (dashed line) and the 
total quark mass m (solid line) at /j = 0, 77 = 0.3 and iV c = 3. 



ing. In the low temperature domain T < 0.4 T c , there is 
a good linear relation between m pg and T. With increas- 
ing temperature, while the order parameter is reduced, 
the pseudogap goes up and dominates the quark mass 
near the transition temperature. Due to the large pseu- 
dogap around the transition, the effective quark mass at 
T c = 146 MeV is about 80% of its value in the vacuum. 
This is in contrary to the mean field picture where the 
quark mass is equivalently considered as the order pa- 
rameter of chiral phase transition and then keeps zero in 
the chiral symmetry restoration phase. 
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FIG. 3: The scaled quark mass m/T as a function of scaled 
temperature T/T c at /1 = 0, r\ = 0.3 and N c — 3 in the region 
above but close to the critical temperature T c . 

Above the transition temperature T c , the order param- 
eter m sc vanishes, but the pseudogap m pg is generally 
not zero in the temperature region close to T c . As we 
mentioned above, the approximation (|19p is in principle 
not valid above the transition temperature. However, for 
temperatures not much larger than T c , it can be used to 
predict the quantitative behavior of pair fluctuations. In 
Fig|3l we calculate the quark mass m in the tempera- 
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FIG. 4: The phase diagram in the T — fj, plane with 77 = 0.3 
and N c — 3. The solid and dashed lines are, respectively, the 
phase transition lines with and without considering pair fluc- 
tuations. The black dots indicate the corresponding tricritical 
points. 



ture domain T c < T < T c mf = 1.3T C , where T c mf = 190 
MeV is the critical temperature at [i = in mean field 
approximation. In this temperature region, the quark 
mass comes purely from the pseudogap. For tempera- 
tures T > 1.2T C , there is approximately a linear relation 
between the quark mass and temperature, m oc T. This 
pseudogap induced quark mass can explain the fact why 
any thermodynamic function can not reach its Stefan- 
Boltzmann limit found by lattice calculations 

We now turn to the case of finite chemical potential. 
The quark chemical potential leads to a mismatch be- 
tween the quark and antiquark Fermi surfaces, like the 
magnetic field which results in a Zeeman splitting be- 
tween the spin-up and spin-down electrons in a supercon- 
ductor. The phase diagram in the T — fx plane is shown 



in Figf4] where T c = 146 MeV is the critical temperature 
at fi = and m(0) = 312 MeV is the quark mass in 
the vacuum. The critical temperature is reduced by pair 
fluctuations, but the pseudogap effect on the phase tran- 
sition decreases with increasing chemical potential. At 
zero temperature, the two phase transition lines with and 
without considering pair fluctuations coincide. In both 
cases, the phase transition is of second order at low chem- 
ical potential and first order at high chemical potentia. 
The tricritical point is shifted from (T, /1) = (78, 285) 
MeV in mean field approximation to (45, 300) MeV when 
pair fluctuations are taken into account. In between the 
two phase transition lines, there exist quarks and pairs, 
like a superconductor in a magnetic field [l8j . 

It is necessary to note that, the pairs considered in 
the GqG scheme are not the real collective excitation 
modes [H}, they just reflect the beyond-BCS effect on 
the fermion propagator. The real collective modes, es- 
pecially the Goldstone modes, are constructed by the 
dressed fermion propagators Gl,r [IH . 

In summary, based on the analogy between the mecha- 
nisms of chiral symmetry breaking and superconducting, 
we have extended the BCS-BEC crossover theory in the 
NJL model to including qq pair fluctuations at finite tem- 
perature and density. By using the BEC condition in the 
symmetry breaking phase, the fluctuations are treated 
as a quark pseudogap which reduces the critical temper- 
ature of chiral phase transition. While the NJL model 
lacks confinement and gluon degrees of freedom, the qq 
pair fluctuations induce a large pseudogap above the crit- 
ical temperature which can be used to explain the large 
thermal quark mass observed in lattice QCD. 
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